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a 
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flux  density 

Pressure 

Electromagnetic  component  of  total  pressure 
Pressure  1n  separator  tank 
Viscous  component  of  total  pressure 


v111 


I 


Electrical  output  power 
Mechanical  Input  power 
Ohmic  power  loss 
Viscous  power  loss 
Volumetric  flow  rate 
Radius 

Inner,  outer  electrode  radii 

Normalized  radius  ratio,  r/r^ 

Internal  resistance  of  equivalent  circuit 

Load  resistance 

Magnetic  Reynolds  number,  p^jCtuL 

Velocity 
Velocity  of  gas 
Velocity  of  1 Iquid 
Load  voltage 

Open-circuit  voltage  of  equivalent  circuit 

Void  fraction 

Generator  efficiency 

Pump  efficiency 

Absolute  viscosity 

Permeability  of  free  space 

Density 

Electrical  conductivity 
Mixture  quality,  mg/Orig  + m^^) 

Cylindrical  coordinate  system 


lx 


DEFINITIONS 


NAME  SYMBOL 

Bubble  Parameter  C' 

d 

Load  Factor  F 

Void  Fraction  a 

Mixture  Quality  x 

Slip  Ratio  K 


DEFINITION 


Dnpirical  ‘drag  coefficient,'  based  on 
experimental  data,  used  to  calculate 
slip  force  between  gas  and  liquid. 

Ratio  of  load  (electrode)  voltage 
to  Induced  voltage,  V|^/Uj^aB. 

Ratio  of  volume  occupied  by  gas  to 
total  volun.e,  or  flow  cross-sectional 
area  occupied  by  gas  to  total  cross- 
sectional  area. 

Ratio  of  gas  mass  flow  rate  to  total 

mass  flow  rate,  m„/(m„  + mj. 

9 9 t 

Ratio  of  gas  to  liquid  velocity. 
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SUMMARY 


The  revised  ambient-temoerature  NaK-nltrogen  facility  is  described.  The 
maximum  liquid  flow  rate  and  generator  inlet  pressure  are  10.9  kg/s  (200  gpm) 
and  1.48  MPa  absolute  (200  psig),  respectively,  compared  with  the  previous 
values  of  6 kg/s  (110  gpm)  and  0.72  MPa  absolute  (90  psig).  Satisfactory  loop 
operation  has  been  obtained,  and  new  experiments  with  the  second  diverging-channel 
generator  were  completed.  The  principal  experimental  results  were  a higher 
power  density  for  the  same  generator  operating  conditions,  and  an  apparent 
tendency  for  the  efficiency  to  improve  more  with  increasing  quality  at  higher 
velocities  than  at  lower  velocities. 

An  evaluation  of  an  annular  generator  geometry  is  presented.  The  advantages 
and  disadvantages  of  the  geometry  are  described,  the  equations  developed,  and 
solutions  obtained  for  three  cases— constant  velocity  and  no  armature  reaction, 
laminar  flow  with  no  armature  reaction,  and  armature  reaction  with  constant 
velocity.  Numerical  examples  show  that:  (1)  the  attainable  terminal  voltages 
appear  to  be  very  low,  (2)  flow  reversal  and  large  viscous  loss  occur  at  or 
below  the  desired  power  densities,  and  (3)  armature  reaction  effects  are  important 
and  compensation  techniques  appear  impractical.  Thus,  this  annular  geometry  does 
not  appear  attractive  for  either  generator  or  pump  operation. 

The  initial  steps  in  the  program  to  produce  and  evaluate  liquid-metal  foams 
are  described.  The  future  directions  of  the  experimental  generator  program, 
including  foams,  are  discussed. 
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I.  INTRODUCTION 


The  two-phase  1 1qu1d-meta1  MHD  generator  program  combines  experimental  and 
analytical  studies  In  an  effort  to  better  understand  the  ohyslcal  processes 
occurring  In  two-phase  generators  and  to  establish  the  limits  of  generator 
performance.  For  the  oast  six  years  this  program  has  been  sponsored  by  the 
Power  Branch  of  the  Office  of  Naval  Research. 

The  generator  experiments  are  conducted  with  an  ambient-temperature 
NaK-N2  test  facility  consisting  of  a two-phase  mixer,  MHD  generator,  magnet,  and 
required  auxiliary  loop  equipment.  The  original  facility,  used  In  all  tests 
through  July  1975,  has  been  described  elsewhere.^  The  revised  facility,  com- 
pleted this  year.  Is  described  In  Section  II. 1 of  this  report.  Earlier  experi- 
ments on  generators  with  varying  geometries  have  been  performed  and  the  results 

934 

presented  In  previous  annual  reports.  ’ ’ 

Analytical  models  of  the  generator  have  been  developed  to  aid  In  the  basic 
understanding  of  the  flow  characteristics  In  the  generator  and  to  provide  a 
(one-dimensional)  computer  code  which  Is  capable  of  approximately  predicting 
generator  performance  from  the  Independent  Input  parameters  such  as  flow  rate, 
magnetic  field  strength,  and  geometry.  This  computer  code  and  the  semi -empirical 
technique  used  to  characterize  the  bubble  parameter  (churn-turbulent  drag  coefficient) 
used  In  the  determination  of  the  Interfaclal  slip  ratio  have  been  reported  In 
detail. 

Previous  generator  experiments  have  used  three  different  channels.  Fig.  I.l. 

The  first  channel  was  designed  for  a silo  ratio  of  unity  as  the  actual  slip  ratio 

was  not  known,  but  higher  slip  ratios  were  observed  experimentally  and  the  liquid 

2 

velocity  was  not  constant,  decreasing  along  the  generator.  Following  this  a 

2 3 

constant-area  channel  was  tested,  and  the  velocity  Increased  along  the  channel.  ’ 
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The  second  diverging  channel,  the  third  channel  teted,  has  a smaller  angle  than 
the  first  channel,  Fig.  I.l,  to  yield  a close-to-constant  liquid  velocity  at  the 

4 

design  point.  As  described  In  the  previous  annual  report.  Improved  generator 
performance  was  obtained  with  this  channel,  and  the  liquid  velocity  was  almost 
constant. 

Previous  experiments  have  been  limited  to  liquid  velocities  of  about 
12  m/s  or  less,  smaller  than  anticipated  for  practical  generators.  In  addition, 
recent  experimental  research  has  demonstrated  that  the  slip  loss  due  to  the 
velocity  difference  between  the  gas  and  the  liquid  can  be  decreased  by  increasing 

3 

the  liquid  velocity.  For  these  reasons  the  facility  has  been  revised  to  allow 
higher  liquid  velocities  or  flow  rates  and  higher  pressures.  In  this  report  the 
revised  facility  Is  described  along  with  the  Initial  data  obtained  with  the  second 
diverging  channel  In  the  revised  facility.  Including  some  high-velocity  runs. 

The  test  program,  described  In  Section  II,  consists  of  the  facility  revisions, 
the  NaK-nItrogen  tests,  and  the  Initial  foam  experiments.  The  results  of  the 
generator  tests,  presented  In  detail  In  Section  III,  Include  a comparison  of 
runs  using  the  same  generator  with  the  original  and  with  the  revised  facilities. 

An  evaluation  of  an  annular  generator  geometry  Is  In  Section  IV.  The  conclusions 
and  suggestions  for  future  work  are  given  In  Section  V. 
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II.  TEST  PROGRAM 

The  experimental  program  Is  summarized  1n  this  section.  The  revised  facility, 
designed  for  generator  tests  at  higher  liquid  velocities  and  pressures.  Is  described 
In  Section  II. 1.  The  program  for  the  NaK-nItrogen  generator  Is  In  Section  II. 2, 
and  the  results  are  described  In  Section  III.  The  Initial  exoerlments  using 
surface-active  agents  to  produce  a foam  are  summarized  In  Section  II. 3. 

II.l  Revised  Ambient-Temperature  NaK-NItrogen  Facility 

The  maximum  liquid  flow  rate  has  been  Increased  from  6 kg/s  to  10.9  kg/s 
(110  gpm  to  200  gpm),  thus  allowing  generator  tests  at  velocities  anticipated 
to  be  typical  of  practical  generators.  As  Increasing  the  flow  rate  Increases 
the  pressure  difference  across  the  generator,  the  available  generator  Inlet 
pressure  was  also  Increased  from  about  0.72  to  1.48  MPa  absolute  (90  to  200  psig). 

The  original  loop  was  constructed  primarily  of  two-inch  pipe  and  had 
several  two-inch  globe  valves  that  caused  severe  flow  restrictions.  In  addition 
the  pressure  tanks  did  not  have  a high-enough  pressure  rating  for  some  of  the 
proposed  tests.  Thus,  the  loop  was  completely  dismantled  and  rebuilt,  as  shown 
In  Fig.  II.l  for  the  NaK  flow  circuit.  The  major  changes  In  the  NaK  circuit  were: 

1.  Two  new  canned-rotor  NaK  pumps  were  added  In  series  with  the  existing 
pump.  The  three  pumps  produce  a gross  pressure  difference  of  1.57  MPa  (228  psi) 
at  10.9  kg/s  (200  gpm).  For  experiments  where  the  total  pump  capacity  Is  not 
required,  one  or  two  pumps  may  be  bypassed. 

2.  The  existing  seoaratlon  and  supply  tanks  were  reolaced  by  a single, 
larger  tank  with  a rated  pressure  of  1.31  MPa  absolute  (175  psig).  It  has 
Internal  baffles  enabling  It  to  separate  twice  the  NaK  and  nitrogen  flow  rate 
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of  the  old  tanks.  A vortex  suppressor  consisting  of  a three-inch  inner  pipe 
surrounded  by  a five-inch  outer  pipe  so  that  the  NaK  flows  up  the  annulus 
between  the  two  pipes  and  down  the  inner  pipe  into  the  line  to  the  pumps  was 
Installed  in  the  tank. 

3.  The  calibration  tank  was  replaced  by  a larger  one  of  210  kg  NaK 
capacity  and  a rated  pressure  of  2.14  MPa  absolute  (295  psig). 

4.  The  150  psi  NaK  filter  housing  was  replaced  by  a 300  psi  housing. 

i 

5.  The  air-cooled  heat  exchanger  was  replaced  with  two  water-cooled  heat 
exchangers  connected  in  parallel.  This  eliminated  the  pressure  loss  of  the 
previous  heat  exchanger,  and  provides  sufficient  cooling  capacity  for  year-round 

operation.  (In  fact,  in  warm  weather  the  loop  can  now  be  made  cool  enough  to  I 

I 

condense  moisture  on  the  NaK  pipes.  Before,  operation  had  to  be  curtailed  in 
warm  weather.) 

6.  Most  of  the  two-inch  piping  was  eliminated,  and  the  NaK  flow  control 
is  now  provided  by  a three-inch  air-operated  valve  rather  than  a two-inch  manual 
globe  valve.  All  non-control  valves  are  ball  valves,  except  for  a four-inch 
Y-pattern  globe  valve  in  the  inlet  line  to  the  first  pump. 

7.  Most  new  components  are  carbon  steel  as  stainless  steel  is  not 
required  for  a room- temperature  NaK  facility. 

8.  All  kerosene  cooling  was  eliminated,  and  the  magnet  and  load  resistor 
are  cooled  by  water.  This  prevents  the  possibility  of  a combined  NaK  and  kerosene 
fire. 

9.  Some  minor  NaK  carryover  with  the  nitrogen  from  the  separator  tank  | 

was  observed.  A three-inch  commercial  secondary  separator  has  been  ordered  | 

and  will  be  Installed  in  the  six-inch  exhaust  line.  j 
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Changes  completed  tn  the  nitrogen  circuit  include; 

'i 

1.  The  nitrogen  purifier  (NaK  bubbler)  tank  was  hydrostatically  tested 

.i 

to  increase  its  oressure  rating.  ij 

2.  An  electrically-operated  valve  to  shut  off  the  nitrogen  supply  in 
case  of  an  emergency  and  a new  pressure  relief  valve  were  installed. 

3.  Some  pipe  sections  and  valves  were  replaced  to  reduce  the  number  of 
bends  and  fittings,  and  the  pipe  length,  and  thus  reduce  the  pressure  drop. 

4.  A proportional  controller  for  the  nitrogen  exhaust  valves  was  added 
to  hold  the  pressure  constant  in  the  separator  tank. 

Extra  nitrogen  heater  capacity  will  be  required  in  the  future  for  the  higher 
flow  rates. 

The  modified  loop  allows  two  modes  of  operation  that  were  not  possible 
before.  These  are: 

1.  Operation  at  up  to  10.9  kg/s  (200  gpm)  NaK  flow  with  a corresponding 
increase  in  the  nitrogen  flow  to  determine  the  effect  of  liquid  velocity  on 
slip  and  generator  efficiency.  Extra  capacity  is  available  and  operation  at  up 
to  18  kg/s  NaK  flow  is  possible. 

2.  Variable  back  pressure  at  the  generator  exit,  up  to  1.14  MPa  absolute 
(150  psig)  at  reduced  flow  rates,  to  determine  the  effect  of  back  pressure  on 
generator  performance.  In  this  mode  only  one  or  two  pumps  may  be  used  to  avoid 
excessive  pressure.  Operation  at  low  back  pressures,  down  to  about  0.12  MPa 
absolute  (17  psia),  is  also  possible,  and  this  allows  larger  pressure  ratios 
across  the  generator.  (The  old  facility  was  limited  to  pressures  above  about 
0.24  MPa  absolute  due  to  the  high  pressure  drop  between  the  separator  tank  and 

the  suction  side  of  the  pump.)  J 
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Several  safety  features  are  Incorporated  Into  the  modifications.  The 
piping  Is  designed  for  2.0  MPa  absolute  (275  psig)  and  the  separation/storage 
tank  for  1.31  MPa  absolute  (175  psig).  Pressure  switches  stop  or  prevent  from 
starting  one,  two,  or  all  three  pumps  If  pre-set  pressures  In  the  separation 
tank  or  high-pressure  piping  are  exceeded.  The  separation  and  nitrogen  purifier 
tanks  are  protected  by  safety  valves,  as  before,  but  with  higher  pressure  ratings. 

Emergency  switches  to  turn  off  the  NaK  pumps  and  the  nitrogen  supply  were  placed 
on  the  main  control  panel  and  adjacent  to  two  doors.  A future  remove  valve 
for  the  cooling-water  supply  has  been  provided  for  In  the  control  panel. 

11. 2 NaK-NItrogen  Generator  Experiments 
The  significant  parameters  for  all  of  the  experimental  runs  for  the 

period  from  August  1975  through  September  1976  are  summarized  In  Appendix  A. 

The  runs  were  chosen  to  give  operating  experience  with  the  revised  facility, 
determine  the  limits  of  the  facility,  and  allow  comparison  of  measurements  with 
the  original  and  revised  facilities  to  see  If  there  were  any  changes.  The  data, 
as  discussed  In  Section  III,  showed  a significant  Improvement  In  generator 
performance  with  the  revised  facility. 

11. 3 Foams 

The  slip  loss  due  to  the  relative  phase  velocities  Is  presently  one  of  the 
major  loss  mechanisms  In  two-phase  liquid-metal  MHD  generators,  especially  at 
the  high  void  fractions  required  for  efficient  cycle  operation.  One  method  of 
reducing  this  loss  Is  to  produce  a more-stable  foam  flow  by  means  of  surface-active 
agents.  The  potential  gains  In  generator  and  cycle  performances  make  this 
approach  very  attractive.  Thus,  a program  to  evaluate  surface-active  agents 
has  been  Initiated.  The  major  points  are: 

I 
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1.  Locate  suitable  agents  for  the  liquid  metals  of  practical  interest. 

2.  Determine  the  techniques  required  for  adequate  separation  of  the 
phases  after  the  generator. 

3.  Perform  tests  in  an  MHD  generator. 

There  is  a lack  of  data  in  the  literature  characterizing  the  preparation  and 
properties  of  a liquid-metal  foam  of  the  desired  kind,  although  the  available 
literature  on  the  characterization  of  aqueous  and  organic  forms  has  given  some 
insight  into  the  preparation  of  a liquid-metal  foam.  In  this  literature  there 
was  frequent  reference  to  the  use  of  surfactant  additives,  and  surface  tension 
and  viscosity  appeared  to  be  the  relevant  physical  properties  of  interest. 
Therefore,  scoping  studies  focused  on  the  effect  of  surfactants  as  liquid-metal 
foam  inducers  and  the  potential  for  fine  particle  suspensions  to  act  as  foam 
promoters  have  been  initiated.  Viscosity  was  selected  as  the  measured  physical 
property  for  use  in  correlation  of  systems  of  differing  bubble  behavior.  In 
later  studies  measurements  of  surface  tension  will  be  included. 

An  experimental  facility  for  studying  surfactants  in  room-temperature 
NaK  has  been  built.  Initially,  the  emphasis  will  be  on  visual  observations 
in  Pyrex  vessels  of  the  foaming  behavior  when  nitrogen  is  bubbled  through. 

A glovebox  was  opened,  cleaned  thoroughly,  and  a new  balance  and  the  gas  system 
for  handling  ultra-high-purity  nitrogen  were  installed.  Safety  of  the  NaK 
operations  was  reviewed.  An  apparatus  for  transferring  NaK  was  designed  and 
built. 

While  waiting  for  the  high-purity  NaK,  experiments  were  done  with  mercury 
to  obtain  experience  with  the  facility  and  to  develop  a "feel"  for  foam  behavior 
in  liquid  metals,  and  encouraging  results  were  obtained.  The  mercury  was  held 
in  a Pyrex  vessel.  An  upturned  tube  made  of  copper  or  of  polyethylene  having 
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a 0.054  cm  orifice  was  used  as  the  gas  bubbler.  Visual  observations  were  made 
of  the  system  behavior  and  a flow  rate  of  5 cc/mln  was  adopted  as  a standard 
condition.  Differences  In  bubble  size  were  noted  when  the  bubbles  were  formed  from 
the  polyethylene  tube  as  compared  with  those  formed  by  the  copper  bubbler. 

The  size  differences  most  certainly  relate,  at  least  In  oart,  to  differences 
In  the  wetting  characteristics  of  the  two  materials.  This  observation  suggests 
that  the  wetting  behavior  between  the  liquid  metal  and  the  mixer  may  be  an 
Important  design  consideration  In  a LMMHD  generator. 

A system  that  was  found  to  show  some  enhancement  of  bubble  lifetime 
was  one  having  one  wt  % cadmium  In  mercury.  Motion  pictures  were  made  of  the 
bubble  action  C45  frames/ second).  The  most  oustanding  enhancement  of  bubble 
lifetime  occurred  when  cesium  at  a concentration  of  3 wppm  was  used  as  the 
surfactant.  For  this  mixture  (Cs  In  Hg)  up  to  seven  coexisting  bubbles  were  seen, 
with  a group  of  six  showing  lifetimes  up  to  10  seconds  after  the  gas  flow  was 
stopped.  Figures  II. 2 and  II. 3,  excerpted  from  the  movie.  Illustrate  the  difference 
In  bubble  behavior  In  the  absence  and  In  the  presence  of  cesium.  Higher  concentrations 
of  cesium  were  also  effective,  but  not  to  the  same  degree.  These  observations 
demonstrate  the  potential  for  developing  a liquid  metal  foam  system  suitable 
for  LMMHD  application. 

The  Initial  experiments  with  NaK  will  focus  on  the  ability  of  fine 
particle  dispersions  to  promote  bubble  lifetime.  Earlier  work  had  Indicated 
that  carbon  black  dispersed  In  NaK  had  apparently  altered  the  surface  properties. 
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III.  GENERATOR  PERFORMANCE 


The  Initial  runs  with  the  second  divergent  channel  In  the  revised  flow 
facility  were  made  for  the  same  conditions  as  with  the  old  facility  to  provide  a 
check  on  the  operation  of  the  facility.  Some  of  these  results  are  summarized 
In  Section  III.l.  A detailed  comparison  for  one  set  of  new  and  old  runs  for  the 
same  operating  parameters  Is  given  In  Section  III. 2.  The  effects  of  reduced 
back  pressure  and  higher  flow  rates  are  discussed  In  Sections  HI. 3 and  III. 4, 
respectively.  A summary  of  all  of  the  experimental  data  for  the  period  from 
August  1975  to  September  1976  Is  given  In  Appendix  A. 

III.l  Repeat  Runs  at  the  Same  Flow  Rate 

The  liquid  velocity,  slip  ratio,  pressure,  local  load  factor,  and  void 
fraction  versus  distance  along  the  generator  are  plotted  In  Figs.  III.l  through 
III. 5 for  one  set  of  open-circuit  runs.  No  major  change  with  respect  to  the 
appropriate  runs  with  the  old  facility.  In  1975,  was  observed.  The  void  fractions 
are  somewhat  higher  than  before,  and  this  Is  probably  connected  with  the  reduced 
slip  and  higher  bubble  parameter  (drag  coefficient),  as  discussed  In  the  next 
section.  There  Is,  as  before,  an  unusually-high  pressure  gradient  at  the  generator 
midsection,  although  this  time  It  seemed  less  pronounced,  and  a pronounced  drop 
In  the  load  factor,  both  at  the  higher  void  fractions. 

The  same  quantities  are  plotted  In  Figs.  III. 6 to  III. 10  for  a generator 
with  a load  resistance  of  R^  * 0.27  mn.  Again,  the  variations  are  similar  to 
the  results  obtained  In  July  1975.  However,  the  pressure  gradient  for  the 
high-vold-fractlon  runs  appeared  to  be  more  uniform  than  previously.  There  Is 
still  a pronounced  decrease  In  the  local  load  factor  with  Increasing  void 
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VOID  FR.ACTION 
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fraction,  probably  due  to  the  onset  of  entrainment  of  liquid  droplets  by  the 
churning  motion  of  the  gas  In  the  generator.  The  fact  that  the  local  slip  ratio 
Is  less  than  unity  at  some  point  Is  not  significant,  this  falls  within  the 
accuracy  of  the  measurement. 


111,2  Comparison  of  New  and  Old  Runs 

Sets  of  data  with  the  old  loop  C7/7/75)  and  the  revised  loop  Immediately 
after  construction  (8/18/76  and  8/24/76)  are  listed  In  Table  III.l.  The  efficiency 
as  a function  of  quality  was  virtually  unchanged,  with  the  exception  of  the 
zero-quality  runs  (see  below).  However,  as  shown  by  the  data  for  a quality  of 
about  0,0058,  the  power  level  (pressure  difference  and  voltage)  was  substantially 
higher  (about  35%  after  correction  for  unequal  magnetic  fields)  for  the  runs 
Inmedlately  after  construction  as  compared  to  the  runs  with  the  old  loop. 

The  runs  with  the  highest  power  level  (August  1976)  had  a lower  measured 
average  velocity  slip  ratio  Y (ratio  of  gas  to  liquid  velocity).  Implying  that 
a more-uniform  two-phase  mixture  existed  In  the  generator,  the  liquid  velocity 
was  higher  (higher  voltage,  current,  and  pressure  drop),  and  the  exit  void  fraction 
(ratio  of  gas  volume  to  total  volume)  was  larger.  Also  the  fluctuation  In  the 
generator  terminal  voltage  was  reduced,  amounting  to  only  a few  percent  of 
the  average  voltage.  (Compare  this  with  the  previous  runs.^) 

The  liquid  velocities  along  the  generator  calculated  from  the  measured 
void  fractions  and  flow  rates  for  the  three  runs  at  about  0.0058  quality  are 
compared  In  Fig.  III. 11.  Note  that  the  differences  In  either  the  variations  of 
or  the  average  liquid  velocities  are  not  sufficient  to  account  for  the  performance 
differences.  An  Important  difference  Is  In  the  average  load  factor  F,  Table  III.l. 
This  Is  higher  for  the  runs  with  the  revised  facility.  Indicating  a higher 
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Table  III.l 
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"effective"  liquid  velocity-electrical  conductivity  combination  for  the  new  runs. 
This  result  Is  expected  from  the  reduced  fluctuations. 

It  is  suspected  that  the  Improvement  in  generator  performance  for  the 
tests  with  the  revised  loop  was  due  to  changes  In  the  surface  properties  of  the 
NaK.  Putting  the  NaK  Into  the  revised  loop  probably  caused  a considerable  Increase 
In  the  Impurities  In  the  NaK,  primarily  iron  oxide  and  Iron  powder  from  the  carbon 
steel.  These  Impurities  could  have  cuased  some  surface  activity  in  the  NaK, 
and  thus  Increased  the  stability  of  the  two-phase  mixture  and  decreased  the  size 
of  the  gas  bubbles  produced  by  the  mixer. 

The  zero-quality  runs  require  special  comment  because  the  reasons  proposed 
above  should  have  no  effect.  Comparing  the  July  1975  and  August  1976  runs,  the 
terminal  voltage  and  efficiency  were  higher  for  the  former  and  the  pressure  drop 
was  lower.  The  measured  void  fraction  profiles  Indicate  that  flow  separation 
occurred  In  at  least  some  of  the  old  zero-quality  runs.  Including  the  July  1975 
case.  Comparing  expected  behaviors  with  and  without  flow  separation,  for  the 
latter  the  liquid  velocity  should  be  less  at  the  exit  end  of  the  generator,  the 
circulating  currents  should  be  higher,  and  thus  the  pressure  gradient  in  the 
upstream  part  of  the  generator  should  be  higher.  (In  the  pure  liquid  case  the 
liquid  velocity  decreases  due  to  the  diverging  channel,  so  that  the  liquid  is 
pumped  In  the  downstream  part  of  the  generator.)  This  predicted  behavior  difference 
Is  consistent  with  the  data.  The  suspected  higher  surface  activity  could  result 
In  Improved  wetting,  and  thus  less  or  no  flow  separation. 

The  calculated  performance  from  the  generator  prediction  code  Is  a terminal 
voltage  of  0.477  V and  an  efficiency  of  0,425  at  1,2  T,  6,0  kg/s  pure  NaK,  and 
0,25  mn  load  resistance.  The  August  1976  runs  yielded  values  close  to  the 
calculated  values. 
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Ill, 3 Reduced  Back  Pressure 


The  liquid  velocity,  slip  ratio,  pressure,  local  load  factor,  and  void 
fraction  along  the  generator  are  shown  In  Figs,  III. 12  through  III. 16  for  a 
back  pressure  only  slightly  above  atmospheric  pressure  In  the  separator  tank. 

In  the  old  facility  It  was  not  possible  to  make  runs  with  separator  tank 
pressures  below  about  0.24  MPa  absolute  (35  psia)  due  to  the  high  pressure  drop 
between  the  separator  tank  and  the  suction  side  of  the  pump.  With  the  improved 
piping  system,  the  lower  limit  on  separator  tank  pressure  Is  that  it  be  high 
enough  to  push  the  nitrogen  through  the  six-inch  line  to  the  atmosphere. 

A reduced  generator  back  pressure  means  that  higher  volumetric  flow  rates 
and  velocities  are  obtained  in  the  generator  with  the  same  mass  flow  rate. 

Also,  lowering  the  back  pressure  automatically  Increases  the  pressure  ratio 
across  the  generator.  This  In  turn  means  that  the  liquid  velocity  will  be 
increasing  along  this  generator  since  It  was  designed  for  a lower  pressure  ratio, 
and  this  should  result  In  reduced  generator  performance. 

The  local  load  factor  seems  to  show  an  Increasing  entrainment  of  liquid 
In  the  gas  along  the  generator.  Fig.  III. 15.  The  void  fraction  and  slip  ratio 
become  quite  large  at  the  generator  exit.  The  expected  Increase  in  the  liquid 
velocity  along  the  generator  Is  observed  for  all  of  the  non-zero  quality  runs. 

Fig.  III. 12. 

A comparison  of  data  for  normal  and  reduced  generator  back  pressures  in 
Appendix  A or  Figs.  III. 6 and  III. 12  shows  that  the  terminal  voltage  Is 
reduced  for  the  lower-back-pressure  runs  even  though  the  average  liquid  velocity 
Is  higher.  The  generator  efficiencies  for  both  back  pressures  are  shown  in 
Fig.  111.17.  The  efficiency  Is  lower  for  the  lower  back  pressure,  probably 
because  of  the  varying  velocity  along  the  generator  and  thus  Increased  circulating 
currents. 
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The  ability  to  go  to  lower  pressure  at  the  generator  exit  will  be  useful 
In  the  future  In  obtaining  larger  pressure  ratios  across  the  generator,  I.e., 
values  closer  to  those  anticipated  for  a commercial  generator.  Of  course,  to 
obtain  high  generator  efficiency  for  higher  pressure  ratios  a generator  taper 
designed  for  the  specific  pressure  ratio  Is  required. 

III. 4 Higher  Flow  Rates 

The  local  oarameters  for  a mass  flow  rate  of  11.9  kg/s  are  shown  In 
Figs.  III. 18  to  III. 22.  The  runs  were  made  at  a reduced  magnetic  flux  density, 
0.6  T,  so  as  not  to  exceed  the  pressure  rating  of  the  channel.  The  attainable 
quality  was  limited  by  some  NaK  carryover  with  the  exhaust  nitrogen,  and  this 
necessitated  Installation  of  a secondary  separator.  Section  II. 1.  Note  the 
reduced  slip  ratios  and  higher  load  factors  that  resulted  for  the  higher  velocity. 
These  are  Indicative  of  a potential  Improvement  In  generator  performance  with 
higher  velocities. 

Generator  efficiencies  at  nominal  flow  rates  of  6 and  12  kg/s  with  the 
other  parameters  the  same  are  compared  In  Fig.  III. 23.  Note  that  there  appears 
to  be  a tendency  for  the  efficiency  to  Improve  more  with  Increasing  quality  at 
the  higher  velocity  than  at  the  lower  velocity.  (There  Is  some  question  about 
the  highest  quality  run  at  the  higher  flow  rate.  Further  runs  were  not  made 
before  the  end  of  the  period  covered  by  this  report  at  this  quality  due  to  NaK 
carryover. ) 
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Figure  III. 23  Generator  Isentrooic  Efficiency  vs.  Quality  at  Different 
Flow  Rates 
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IV.  ANALYTICAL  STUDIES  OF  ANNULAR  GENERATORS 


IV. 1 Introduction 

Rectangular  MHD  conduction  machine  geometries  for  generators  and  liquid- 
metal  pumps  have  been  studied  almost  exclusively.  Helical  geometries  are 
suitable  for  low-volume-flow-rate  high-pressure  difference  pumps,  but  not  for 
MHD  generators.  Annular  geometries  are  common  for  Induction  pumps  and  generators, 
but  there  are  no  known  references  on  annular  conduction  machines.  Since  the 
annular  geometry  has  several  advantages,  the  following  analytical  study  was 
undertaken  to  evaluate  Its  practicality. 

The  annular  generator  geometry  considered  consists  of  two  concentric, 
homogeneous  cylindrical  walls  of  different  radii  with  the  conducting  fluid  flowing 
In  the  axial  direction  within  the  annulus.  Fig.  IV. 1.  Advantages  of  this  geometry 
are; 

1.  The  cylindrical  walls  are  the  electrodes,  eliminating  the  need  for 
Insulating  walls  and  the  resulting  fabrication  problems.  Without  the  Insulating 
walls  there  are  no  accompanying  viscous  boundary  layer  and  current  shunting 
effects  as  found  In  rectangular  geometries. 

2.  Homogeneous  cylindrical  shells  are  Ideal  pressure  vessels. 

3.  The  circumferentially-directed  magnetic  field  Is  entirely  contained 
within  the  generator,  stray  fields  are  not  a problem,  and  the  field  Is  produced 
by  a superconducting  or  ordinary  torroldal  winding  without  any  magnetic  Iron. 
Disdavantages  of  this  geometry  are; 

1.  The  electrical  terminal  properties  are  Inherently  low  voltage  and  high 
current  because  of  the  relatively-short  electrode  spacing  and  large  cross-sectional 
area  for  current  flow. 


EXCITING  TOROIDAL  WINDING, 
NIo.r  AMPERE-TURNS 


2.  The  electromagnetic  force  density  has  an  Inverse-square  radial  dependence 
which  may  result  In  a very  nonuniform  velocity  profile. 

3.  The  magnetic  field  generated  by  the  Induced  currents  (armature  reaction) 
may  cause  large  variations  In  the  power  density  along  the  machine,  and  there  Is 

no  apparent  method  to  comoensate  for  this  without  destroying  the  first  advantage. 

4.  The  toroidal  magnet  winding  has  to  be  closed  through  the  end  regions 
of  the  annular  duct,  so  that  the  wires  must  cross  the  flow  area. 

The  first  three  disadvantages  are  analyzed  In  the  following  material  by 
means  of  three  simplified  models.  After  establishing  the  relevant  equations 
and  assumotlons,  a model  with  a constant  velocity  and  no  armature  reaction  Is 
considered  to  establish  practical  regions  of  ooeratlon  In  Section  IV. 3.  A 
laminar  velocity  profile  without  armature  reaction.  Section  IV. 4,  Is  used  to 
study  the  effect  of  the  Inverse-square  radial  dependence  of  the  electromagnetic 
force  density.  The  flow  In  a practical  generator  could  be  either  laminar  or 
turbulent,  depending  on  the  relative  1 ami narl zing  effect  of  the  magnetic  field, 
but  the  laminar  solution  yields  Insight  Into  what  might  happen  even  for  turbulent 
flow.  Finally,  the  axial  deoendences  of  the  magnetic  flux  and  current  densities 
for  a constant  fluid  velocity  are  considered  In  Section  IV. 5.  Numerical  results 
obtained  from  an  annular-geometry  computer  code  are  Included  for  all  three 
models.  The  primary  objective  of  this  study  Is  generator  operation,  but  the 
solutions  apply  equally  to  Dump  operation  and  numerical  oump  cases  are  also 
discussed.  Only  pure-liquid  machines  are  considered,  but  the  conclusions  are 
valid  for  two-phase  liquid-metal  and  plasma  machines  as  well.  The  velocity- 
profile  solution  also  has  Important  Implications  for  annular  Induction  pumps 
and  generators. 


-45- 


IV. 2 Equations  and  Assumptions 

All  models  considered  are  qoverned  by  the  same  general  equations  of 
fluid  mechanics  and  electromagnetics.  The  assumptions  made  to  solve  the 
equations  are: 

1.  Steady-state  conditions, 

2.  Incompressible,  Newtonian,  non-polarizable  fluid,  I.e.,  pure 
liquid, 

3.  No  fluid  sources  or  sinks  In  the  generator, 

4.  Negligible  gravitational  effects, 

5.  One-dimensional  axial  flow  with  circular  symmetry  about  the  axis 
of  the  cylindrical  walls,  and 

6.  The  magnetic  permeability  Is  that  of  free  space  everywhere. 

With  these  assumptions,  the  basic  relevant  equations  are 

V X 5 = 'S,  (IV. 1) 

0 

= o(t  + u X t),  (IV. 2) 

V X t = 0,  (IV. 3) 

7 . S = 0,  (IV. 4) 

7 • u = 0,  (IV. 5) 

and  p(u  • 7)u  = -7p  + pV^u  + ^ x (IV. 6) 

where  the  variables  are  defined  In  the  nomenclature. 

In  cylindrical  coordinates  with  B.,  u,,  J„,  and  as  the  only  vector 
components  the  equations  become 


If  “ 0* 


(IV. 7) 


3B. 


(IV. 8) 


(IV 
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(rv.9) 

(IV. 10) 

(IV. 11) 

(IV. 12) 
(IV. 13) 


The  maqnetic  flux  density  B-  1s  most-easi ly  determined  from  the  Integral 

u 

form  of  Eq.  IV. 1 evaluated  over  a concentric  circle  about  the  axis,  or 


(I  = 2Trr  B-  = u ? • n da. 

9 0 


(IV. 14) 


C S 

For  the  applied  field,  the  total  field  without  armature  reaction. 


r»  _ ''^^exc 




(IV. 15) 


where  Is  the  ampere-turns  of  the  exciting  torroldal  winding.  Fig.  IV. 1. 

The  electrical  terminal  properties,  the  electrode  or  load  current  1^ 
and  the  electrode  voltage  difference  or  load  voltage  Vj^,  are  given  by 


(IV. 16) 

(IV. 17) 
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1 


where  L Is  the  length  of  the  machine,  and  and  are  the  radii  of  the  Inner 

1 0 

and  outer  cylinders  (electrodes). 

The  general  equations  for  the  powers  are  collected  here  for  reference. 
The  ohmic  loss  In  the  fluid  Is 


/ ^ 


rdrdedz. 


(IV. 18) 


The  electrical  output  power  with  end  loss  neglected  (negative  for  a pump)  Is 
Pg  = (IV. 19) 

The  mechanical  Input  power  (negative  for  a oump)  Is 

’’m  ■ '’(‘Pm  * ‘Pv)-  ('''.20) 

where  Q Is  the  volumetric  flow  rate  of  the  working  fluid,  Ap^  Is  the  electro- 

60) 

magnetic  pressure  difference. 


fpLu  2 


(IV. 21) 


Is  the  (non-MHD)  viscous  pressure  difference  In  terms  of  the  Darcy  friction  factor 
f,  and  = 2(rp  - r^)  Is  the  hydraulic  diameter.  The  efficiencies  for  generator 
and  pump  operation  are  defined  as 


® - "I 

- and  n.,  = K- 


e 


(IV. 22) 


It  Is  convenient  to  use  the  normalized  or  non-dimensional  radius 


R = r/r^ 

In  the  analysis.  At  the  outer  electrode. 


(IV. 23) 


(IV. 23a) 
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IV. 3 Constant  Velocity  with  no  Armature  Reaction 


Constant  fluid  velocity  with  no  amature  reaction  (the  ideal  case) 
is  assumed  to  establish  the  parameters  required  for  practical  applications.  The 
required  equations  are  determined,  and  numerical  examples  calculated  for  both 
generator  and  pump  operation. 

IV. 3.1  Analysis 

The  load  or  terminal  current  is 


1l  * -2nrU^  (IV. 24) 

from  Eq.  IV. 16.  The  terminal  voltage  is  calculated  from  Eq.  IV.17,  where 
E^  is  obtained  by  substituting  Ij^  and  from  Eq.  IV. 15  into  Eq.  IV. 9 and 
solving  for  E^.  The  result  is 


“z  ".'"exc 


(IV. 25) 


in  terms  of  the  normalized  radius.  This  can  be  interpreted  as  the  equivalent 


circuit  of  Fig.  IV. 2,  which  has  the  equation 

'l  - 'oc  - 'l'i- 
The  open-circuit  voltage  is 


{IV. 26) 


“z  '‘o'^'exc  l"«o 


(IV. 27) 


and  the  internal  resistance  is 


Ri  = lnRQ/(27rLo).  (IV. 28) 

In  terms  of  the  magnetic  flux  density  at  the  inner  radius,  B^,  V^^  can  be 


written  as 


V_  » u^B.r^lnR  . 
oc  z 1 1 0 


{IV. 27a) 
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The  ohmic  power  loss  in  the  fluid,  from  Eq.  IV. 18,  is 

(IV. 29) 

The  electrical  output  power,  Eq.  IV. 19,  is  written  in  terms  of  V^  as 
Pg  • oLVL(y  u^  Nlg^^  - 2n  V^/lnR^j)  (IV.30) 

using  Eq.  IV. 25  or  Eq.  IV. 26  to  eliminate  I|^.  In  terms  of  the  load  factor,  F, 
defined  as 


^ ■ vv- 

Eqs.  IV. 29  and  IV.30  can  be  written  as 


and 


(IV. 31) 

(IV. 29a) 


P.  ■ (Vjc“/R,)  F(1  - F).  (IV.30.) 

The  calculation  of  is  more  complicated  because  the  electromagnetic  pressure 
gradient  is  not  constant  over  the  channel  cross  section.  To  calculate  Ap^^^, 
the  Lorentz  force  density  ^ x ^ Is  integrated  over  the  total  channel  volume, 
yielding 


■ -“N'exc 

The  minus  sign  means  that  F^  acts  in  the  -z  direction  or  opposes  the  fluid  motion. 

cfn 

Then, 


AP 


em 


w NI  I,  InR^ 

0 exc  L 0 

2,r2  r^2  (Rq2  . 1) 


(IV. 33) 


is  the  average  difference.  The  mechanical  input  power  is  calculated  from 
Eqs.  IV. 20,  IV. 21,  and  IV. 32,  and  the  efficiency  from  Eq.  IV. 22.  If  the 
viscous  pressure  difference  is  negligible,  then 

Pm  - 
and 

Og  - F (IV. 35) 

as  in  rectangular  MHO  machines. 


i 
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IV. 3. 2 Numerical  Examples 

The  most  significant  parameter  Is  the  obtainable  terminal  voltage,  as 

Vj^  must  be  large  enough  for  either  efficient  dc  to  ac  Inversion  or  homopolar 

motors.  Since  the  load  factor  F must  be  close  to  unity  for  efficient  generator 

or  Dump  operation.  Eg.  IV. 35,  V^^  Is  given  In  Table  IV. 1 for  values  of  r^,  R^, 

and  u,.  Note  that  R..  Is  limited  to  be  much  less  than  1.1  at  reasonable  power 
z 0 

densities  to  prevent  the  possibility  of  flow  reversal,  as  shown  In  Section  IV. 4. 
The  velocity  Is  limited  to  around  25  to  50  m/s  by  viscous  loss.  None  of  the 
cases  shown  yields  adequate  voltage,  even  at  = 6 T and  r^  = 5 m.  Higher 
voltages  are  obtainable  with  several  generators  connected  electrically  In  series, 
except  that  the  output  power  per  generator  (see  Table  IV. 2)  may  be  too  large 
to  use  several  generators  In  one  power  plant. 

The  powers  are  shown  In  Tables  IV. 2 as  a function  of  the  load  factor  F 
for  one  of  the  cases  from  Table  IV. 1.  For  an  acceptable  efficiency  with  all 
losses  Included  (end  losses  are  Ignored  here)  F will  most  likely  be  between 
0.9  and  0.95,  leading  to  a generator  power  on  the  order  of  500  MVI  per  meter 
of  length  and  V^^  around  7 V at  r^  = 1 m and  = 3 T.  Since  the  power  scales 
as  Vq^2  (o,.  |r.^2  8^2)  3^  constant  R^  and  F,  Eq.  IV. 30a,  Increasing  V^^ 

results  In  Increased  and  no  possibility  of  a series  connection  to  raise  the 
useable  voltage.  Finding  a good  compromise  between  generator  size,  power  level, 
and  terminal  voltage  Is  very  difficult  for  a rectangular  generator,  and  still 
more  d1ff1cult,1f  not  Impossible,  for  large  annular  generators. 

There  are  no  applications  for  electromagnetic  pumps  of  the  power  levels 
shown  In  Table  IV. 2.  Thus,  acceptable  voltage  levels  may  be  unattainable  for 
annular  conduction  pumps  since  the  size  and  voltage  must  be  reduced  from  the 
values  of  Table  IV. 2.  The  low  voltages  shown  Table  IV. 2 are  not  surprising, 
they  result  from  the  relatively-short  electrode  spacing  mentioned  In  Section  IV. 1 
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Table  IV. 1 


Open-Circuit  Voltage  and  Excitation  Ampere-Turns  as  a Function  of 
Electrode  Radii,  Magnetic  Field  Strength  at  r^,  and  Constant  Velocity 


"1 

[m] 

>^0 

“i 

[T] 

“z 

[m/s] 

Nlexc 
[10^  AT] 

V 

Cv] 

1 

1.1 

3 

25 

1.5 

7.15 

1 

1.1 

3 

50 

1.5 

14.3 

1 

1.1 

6 

25 

3.0 

14.3 

1 

1.2 

3 

25 

1.3 

13.7 

3 

1.1 

3 

25 

4.5 

21.4 

3 

1.01 

3 

25 

4.5 

2.24 

3 

1.1 

6 

25 

9.0 

42.9 

3 

1.2 

3 

25 

4.5 

27.3 

5 

1.1 

3 

25 

7.5 

35.7 

5 

1.01 

3 

25 

7.5 

3.73 

5 

1.1 

6 

25 

15.0 

71.5 

5 

1.1 

3 

50 

7.5 

71.5 

5 

1.2 

3 

25 

7.5 

68.4 

Table  IV. 2 


Powers  per  Meter  Length  for  r^ 


Fluid  = NaK 
B.  = 3 T 
Uj  = 25  m/s 
= 1.5  X 
7.15  V 


V 

oc 


lO'  A- 


1 m,  = 1.1,  and  Constant  Velocity 

p = 867  kg/m^ 
a « 2.62  X 10°  mhos/m 

°h  = 

f*  0.022 


F 

''l 

Pe 

Pm 

Useful  Power 
Density 
[MW/m3] 

"g 

[V] 

[MA] 

[MW] 

[MW] 

[X] 

Gen. 

0.8 

5.72 

247 

1410 

1770 

2140 

79.9 

Gen. 

0.95 

6.79 

61.8 

420 

442 

636 

94.9 

Pump 

1.1 

7.86 

-123 

-971 

-882 

1340 

90.9 

Pump 

1.3 

9.29 

-370 

-3440 

-2650 

4010 

76.9 

F 


i 

V 

1 

1 

I 


( 

S 

r 

! 

I 

} 

( 

f 


IV. 4 Fully-Developed  Laminar  Flow  with  no  Armature  Reaction 

In  this  section  the  effect  of  the  radial  dependence  of  the  electromagnetic 
force  density  Is  analyzed.  Fully-developed  laminar  flow  Is  treated  because 
analytical  solutions  can  be  obtained.  The  laminar-flow  solution  provides  a 
guide  for  the  behavior  of  turbulent  flows. 


I 


IV. 4.1  Analysis 

First  the  velocity  profile  u^(r)  must  be  determined.  The  electromagnetic 
force  density  0 B.  can  be  written  as 

r 0 


Jr«9 


% ^^exc 
(27rr)2L 


•*"1  ^1 
2irr2L 


(IV. 36) 


from  Eqs.  IV. 15  and  IV. 24.  The  pressure  gradient  can  be  expressed  as  3p/3z  = -Ap/L, 
where  Ap  equals  the  Inlet  minus  the  exit  pressures.  Note  that  Ap  Is  Independent 
of  r and  z for  fully-developed  flow  and  Is  positive  for  generator  operation. 
Substituting  Into  Eq.  IV. 13  and  solving  for  u^{r}  with  the  boundary  conditions 
that  u^  = 0 at  the  walls  yields 


u^CR)  * [-In  R InCR^/R)] 


'l 

4it  Ly 


^f^lnR-  (R^-l) 


r^2  AP 

"irnr 


(IV. 37) 

In  terms  of  the  normalized  radius  R.  Note  that  this  Is  not  a Hartmann  profile, 
and  there  Is  no  parameter  like  the  Hartmann  number  governing  the  shape  of  the 
velocity  profile. 

Next  Vj^  Is  calculated  using  the  same  procedure  as  In  Section  IV. 3.1 
except  that  u^(R)  Is  given  by  Eq.  IV. 37.  The  result.  In  terms  of  the  equivalent 
circuit  of  Fig.  IV. 2 or  Eq.  IV. 26,  Is 
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(IV. 38) 


and 


V • R,  - 1)  + 1 * In  y 


In  R.  r,2  B,2(ln  R.)’ 
0 . 1 1 0 

2it  La  24TrLu 


Ap 

m 


(IV. 39) 


Note  that  1s  not  the  physical  Internal  resistance  In  Rjj/2itLct.  The  second 
term,  which  Is  Independent  of  o,  comes  from  the  specified  constant  pressure 
difference.  With  Ap  fixed,  the  average  fluid  velocity  or  the  volume  flow  rate 
Q decreases  as  1^^  Increases  because  of  the  Increased  (adverse)  electromagnetic 
pressure  gradient.  In  Section  IV. 3.1,  on  the  other  hand,  u^  (or  Q)  was  specified 
and  Ap  had  to  Increase  with  Increasing  Ij^. 

The  ohmic  power  loss  Is  still  given  by  Eq.  IV. 29  since  the  second  term 
of  Eq.  IV. 39  does  not  represent  a physical  resistance.  The  electrical  output 
power  Pg  = Vj^Ij^  Is  expressed  as  a function  of  Ap  and  1^^  using  Eqs.  IV. 26,  IV. 38, 
and  IV. 39.  The  mechanical  Input  power  Is  calculated  from  = QAp,  where  the 
volumetric  flow  rate  Is 


Q 


2Trru^dr. 


Integrating  Eq.  IV. 40  using  Eq.  IV. 37  results  In 


(IV. 40) 


Pm  = tR.^dnR,  - 1)  + 1 + InRg] 


Po^  - 1 

Tnr~ 


r^3  1^  Ap 

5U^ 

irri**  (Rg2  - 1)  Ap2 

sm 


(IV. 41) 


The  viscous  power  dissipation  can  be  calculated  from 
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% - 27rL 


(IV. 42) 


0 3U_ 

u(  3^  rdr 


i 


using  Eq.  IV. 37,  or  from  the  other  powers  and  conservation  of  energy. 


IV. 4. 2 Numerical  Examples 

The  velocity  profiles  for  r^  = 3 m,  * 1 .01 , = 3 T,  V|^  =•  4 V,  NaK, 

and  three  pressure  gradients  are  shown  in  Fig.  IV. 3.  The  powers  for  these 
parameters  and  six  pressure  gradients,  including  the  three  plotted  in  Fig.  IV. 3, 
are  given  in  Table  IV. 3.  Increasing  pressure  gradient  corresponds  to  decreasing 
load  factor  F,  as  shown  in  Table  IV. 3.  Note  that  for  small  pressure  gradients 
(F  close  to  unity)  the  profile  resembles  the  non-electromagnetic  annular 
profile,  the  viscous  and  ohmic  losses  are  small,  and  the  power  density  is  low. 

At  a power  density  of  about  30  MW/m^  flow  reversal  occurs  (Fig.  IV. 3)  and  the 
viscous  power  loss  is  substantial.  Although  the  efficiency  is  very  high,  the 
load  factor  is  so  close  to  unity  that  the  power  density  is  relatively  low  for 
practical  MHO  generators  and  the  end  losses  (not  included)  would  be  serious. 

A similar  behavior  is  observed  for  R « 1.1,  Fig.  IV. 4,  except  that  the  : 

i 

flow  reversal  occurs  at  a lower  pressure  gradient  or  power  density  because  of  i 

the  larger  variation  in  the  electromagnetic  pressure  gradient  across  the  channel  | 

cross  section.  I 

Note  that  the  behavior  of  the  annular  generator  is  different  from  the  | 

I 

convertional  rectangular  generator  in  that  flow  reversal  occurs  and  the  viscous  ! 

j 

losses  increase  sharply  with  generator  loading  (decreasing  F).  Flow  reversal  | 

occurs  for  all  cases  except  very  low  power  densities  so  that  the  viscous  j 

contribution  to  the  total  pressure  gradient  can  balance  out  the  radial  variation  | 

of  the  electromagnetic  pressure  gradient.  In  the  rectangular  generator  | 
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Table  IV. 3 

Powers  per  Meter  Length  and  Efficiencies  as  a Function  of  the 
Pressure  Gradient  for  r^  - 3.0  m,  R.  =•  1.01, 

Bi  - 3 T,  Vl  - 4 V,  and  NaK 


Ap/L 

[N/m^] 

F 

^1 

[*3 

Useful  Power 
Density 
[MM/ni3] 

% 

[MW] 

Pm 

m 

[MW] 

% 

[MW] 

Pr 

[MW] 

6.9  X 10^ 

1.0000 

27.0 

0.00832 

0.00473 

0.0175 

0.0128 

8.44  X 10"’° 

6.9  X 10^ 

1.0000 

92.7 

0.286 

0.162 

0.175 

0.0128 

9.96  X 10"^ 

6.9  X 10^ 

0.9999 

99.2 

3.06 

1.74 

1.75 

0.0144 

1.14  X 10'^ 

6.9  X 10^ 

0.9993 

99.0 

30.8 

17.5 

17.7 

0.171 

1.16  X 10"^ 

6.9  X 10® 

0.9934 

91.2 

308. 

175. 

192. 

15.8 

1.16 

6.9  X 10^ 

0.9379 

50.8 

3080. 

1750. 

3450. 

1580. 

116. 

r ^ 

i 

I 

i 

I 

I 

i 

r 

viscous  power  loss  1s  almost  constant,  except  for  the  slight  Increase  In  average 
velocity  needed  to  hold  the  load  voltage  constant  as  F Is  decreased.  In  the 
cases  considered  for  the  annular  generator  the  viscous  loss  Is  always  much  larger 
than  the  ohmic  loss,  and  the  efficiency  decrease  at  higher  power  density 
(lower  F)  Is  caused  by  the  sharp  Increase  In  the  viscous  loss  with  the  flow 
reversal  and  the  resulting  large  shear  forces. 

The  velocity  profiles  for  two  sets  of  pump  cases,  = 1.01  and  1.1,  are 
shown  In  Figs.  IV. 5 and  IV. 6,  respectively.  The  corresponding  powers  and 
efficiencies  are  given  In  Tables  IV. 4 and  IV. 5.  In  both  cases  flow  reversal 
occurs,  although  at  a much  higher  pressure  gradient  or  power  density  for  the 
Rp  = 1.01  case  because  of  the  smaller  change  In  the  electromagnetic  force  density 
across  the  channel.  Note  that  for  even  higher  pressure  gradients  or  power 
densities  the  devices  become  dissipators,  absorbing  both  electrical  and  mechanical 
Input  powers  and  dissipating  them  as  ohmic  and  viscous  losses.  (A  similar 
effect  will  occur  with  any  pump.  A pump  has  a given  pressure  difference  with 
zero  flow;  If  a higher  pressure  difference  Is  Imposed,  then  reversal  of  the  net 
flow  will  occur. ) 

the  voltage  levels  for  the  pump  cases  are  very  low  because  of  the  short 
electrode  spacing.  Both  cases  considered  have  the  same  average  electric-field 
strength.  The  low  voltages  mean  that  the  device  Is  probably  not  practical. 

The  results  obtained  here  are  only  valid  for  laminar  flows.  However, 
similar  effects  are  expected  with  turbulent  flows.  Including  flow  reversal 
at  high  power  densities. 
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VELOCITY  (m/s) 


1 


NORMALIZED  RADIUS  R 


Fiqure  IV, 5. 


Velocity  versus  Radius  as  a Function  of  the  Pressure  Gradient 
for  a Pump,  ■ 1,01 
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VELOCITY  (m/s) 


Table  IV. 4 


Powers  per  Meter  Length  and  Efficiencies  as  a Function  of  the 
Pressure  Gradient  for  a Pump,  r<  ■ 0.5  m,  R-  « 1,01, 

Bi  > 2 T,  “ 0.05  V,  and  NaK 


Ap/L 

CN/ln’] 

F 

'’P 

[«] 

Useful  Power  -P. 
Density 

[kW/m^J  [KH] 

-Pm 

[kH] 

[ua 

Pr 

m 

-6.9  X 10^ 

1.000 

25 

3.466  0,2166 

0,05473 

0.1619 

1.13  X 10-® 

-6.9  X 10^ 

1.000 

77 

34.66  0,7092 

0.5472 

0.1618 

1.22  X 10-^ 

-6.9  X 10^ 

1.001 

97 

346.1  5.635 

5.465 

0.1626 

7.6  X 10-^ 

-6.9  X 10^ 

1.013 

98 

3,413  54.89 

53.89 

0.2798 

0.7286 

-6.9  X 10® 

1.153 

85 

29,300  547.5 

462.7 

12.34 

72.47 

-6.9  X 10^ 

- 

- 

5,473 

-2,992 

1,222 

7,243 

Table  IV. 5 

Powers  per  Meter  Length  and  Efficiencies  as  a Function  of  the 
Pressure  Gradient  for  a Pump,  r-i  » 0.5,  Rq  ■ 1.1, 

Bi  - 2 T,  Vl  =■  0.5  V,  and  NaK 

AP/L 

[N/m^] 

F 

^P 

L%1 

Useful  Power  -P_ 

Density 

[kW/m3]  [kW] 

-Pm 

[kW] 

Pv 

[kW] 

Pr 

[kva 

-6.9  X 10^ 

1.000 

95 

3.73  0.6151 

0.5849 

0,0301 

8.7  X 10-® 

-6.9  X 10^ 

1.000 

80 

36.30  5.985 

4.796 

1.189 

8.3  X 10-^ 

-6.9  X 10^ 

1.001 

- 

59.68 

-57.44 

117.0 

0.0825 
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IV. 5 Armature  Reaction  with  Constant  Ve1oc'tt\ 

The  effect  of  the  axial  variation  of  the  magnetic  flux  density  and  the 
(electric)  current  density  caused  by  the  current  in  the  fluid  is  considered 
here.  The  three  different  connections  to  complete  the  load  circuit,  shown  in 
Fig.  IV. 7,  are: 

1.  At  the  downstream  end, 

2.  At  the  uDstream  end,  and 

3.  At  both  ends. 

The  three  connections,  their  different  effects  on  the  magnetic  field,  and  the 
resulting  performances  are  discussed  in  the  following  material. 

IV. 5.1  Downstream  Connection 

If  the  load  is  connected  to  the  downstream  ends  of  the  electrodes,  then 
the  z-directed  load  current  flowing  along  the  inner  electrode  is  in  the  same 
direction  as  the  exciting  current  for  a generator.  Fig.  IV. 7.  Thus,  the 
induced  magnetic  field  due  to  the  load  current  reinforces  the  applied  field, 
and  the  total  field  is  every  place  greater  than  the  aonlied  field  (except  at  the 
generator  entrance  where  the  induced  field  is  zero). 

The  apolied  field  is  riven  by  Eq.  IV. 15,  and  the  total  field  is 

NI  ^ 

Bg(r,z)  » J/r,z')dz'  (IV.43) 

0 


from  Eq.  IV. 14  and  conservation  of  charge.  The  electric  field  E^  must  vary 
as  1/r  from  conservation  of  charge  or  Eq.  IV. 8, and  Eq.  IV. 9;  and  in  terms  of  V, , 


Figure  IV. 7 Load  Connections  for  Armature  Reaction  Calculations 


-66- 


from  Eq.  IV. 17.  Substituting  Into  Ohms  law,  Eq.  IV. 9,  gives 


o V. 


Rm 


J^(r.z)  - ^ 


Nlexc 


‘I 


J^(r,z')d2' 


(IV. 45) 


where 


R_  = y o u_  L 
m 0 z 


(IV. 46) 


Is  the  magnetic  Reynolds  number.  Solving  Eq.  IV. 45  yields 


Jy.(r.z) 


"mNW 

~5iirU 


oVl 

rlrTR’. 


0 J 


(IV. 47) 


and,  from  Eqs.  IV. 43  and  IV. 16, 


B9(r,z)  = 


>1  z/L 

urWt-  ^ ^^^.48) 


and 


I 


L 


N'exc 


ZttV, 


0 Z 0 


1). 


(IV. 49) 


From  comparing  Eq.  IV. 49  with  the  equivalent  circuit.  Fig.  IV. 2,  V Is  given 
by  Eq.  IV. 27  for  the  Ideal  case,  constant  velocity  with  no  armature  reaction, 
since  there  Is  no  armature  reaction  for  open  circuit.  The  Internal  resistance 
for  this  case. 


R 


1 


2nLa(e 


(IV. 50) 


Is  always  less  than  the  value  for  the  Ideal  case,  Eq.  IV. 28.  The  reduced  value 


for  R^  occurs  because  the  total  magnetic  flux  In  the  generator  Increases  with 


loading,  and  this  acts  to  compensate  for  the  Increased  ohmic  voltage 
drop. 


r 


t 


1 


I 

I 


The  Lorentz  force  Integrated  over  the  total  channel  volume  Is 


2 - 1)  + (1  - F)(e  -1)3. 


(IV. 51) 

From  this  plus  the  viscous  oressure  drop,  Eq.  IV. 21,  the  mechanical  Input  power 
can  be  calculated.  The  ohmic  power  loss  In  the  fluid  Is 


"r  ■ 1!^  (V  - " - ')■ 


2R_ 


(IV. 52) 


The  remaining  powers  and  efficiencies  can  be  calculated  from  the  equations 
In  Section  IV. 2. 


IV. 5.2  Upstream  Connection 

If  the  load  Is  connected  to  the  upstream  ends  of  the  electrodes,  then  the 
z-directed  load  current  flowing  along  the  Inner  electrode  Is  In  the  opposite 
direction  to  the  exciting  current.  Fig.  IV. 7.  For  this  case  the  induced  magnetic 
field  opposes  the  applied  field,  the  total  field  Is  everyplace  less  than  the 
applied  field  (except  at  the  generator  exit  where  the  Induced  field  Is  zero), 
and  the  reduction  In  the  total  magnetic  flux  In  the  generator  Is  reflected  In  a 
larger  Internal  resistance. 

The  procedure  used  to  determine  the  equations  for  the  upstream  connection 
Is  the  same  as  for  the  downstream  connection.  The  results  are; 


J 
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r 


B,(r.z) 


"'exe 
“ 


I 


J^(r,z')dz'. 


(IV. 43a) 


J^(r,z) 


'^"W  “''l 

2.rL  ■fTiTT  * 

0 . 


(IV. 47a) 


, V V,  R«(z/L-1) 

Z 0 


ZirV. 


NI 


exc  ■ il  u.  In 

0 Z 0 


(1  - 


M ^ 

2irLa(l  - e ) 


(IV. 48a) 
(IV. 49a) 


(IV. 50a) 


and 


■w. 


.-"m, 


"of'exc*^'’  - ■ « ) * 2F(1  - e ")]. 

(IV.SIa) 

Note  that  V^^  1s  unchanged,  and  Is  always  greater  than  the  value  for  the 
Ideal  case,  Eq.  IV. 28. 


IV. 5. 3 Both  Ends  Connection 

With  the  load  connected  to  both  ends  of  the  generator.  Fig.  IV. 7,  the 
generator  Is  split  Into  two  regions  by  the  location  z * L'  where  the  component 
of  the  load  current  flowing  along  the  axis  Is  zero.  The  region  for  z < L' 
behaves  like  a generator  with  the  upstream  connection,  the  region  with  z > L' 
behaves  like  a generator  with  the  downstream  connection,  and  the  total  properties 
can  be  determined  from  the  sum  of  the  two  regions  using  the  equations  from 
Sections  IV. 5.1  and  IV. 5. 2.  For  example,  the  load  current  Is  given  by 
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^LU  ^ ^LD 


NI 


exc 


2ttVl 

p u.  In 
0 z 0 


(e  " - e 


V''!- 

" ). 

(IV. 49b) 


where  Is  based  on  the  total  length  L. 

The  location  of  L'  would  normally  be  determined  by  the  voltage  drop 
along  the  Inner  electrode.  Since  the  electrode  Is  assumed  to  be  an  equipotentlal 
here.  It  Is  not  possible  to  calculate  a value  for  L'. 


IV. 5.4  Numerical  Example 

The  powers  and  efficiencies  for  one  set  of  parameters  are  shown  1n 

I 

Table  IV. 6 for  the  three  connections,  and  the  magnetic  flux  densities  and 
current  densities  are  shown  In  Figs.  IV.3  and  IV. 9.  For  comparison,  the  values 
with  the  armature  reaction  Ignored  are  also  shown.  The  value  = 4.939 

Is  not  large,  yet  the  effects  are  dramatic.  In  fact,  the  downstream  connection 
results  are  not  realistic,  the  pressure  gradient  and  magnetic  flux  density 
both  become  unrealistically  large.  (No  Iron  Is  used,  so  that  magnetic  saturation 
1s  not  a factor.)  The  upstream  connection  Is  realizable,  but  the  power 
density  Is  low.  The  choices  u^  » 15  m/s  and  a ■ 2.62  x 10^  mhos/m  (0.1  times 
a of  NaK,  characteristic  of  a high-vold-fractlon  flow)  were  made  to  restrict 
to  small  values.  The  excitation  was  chosen  to  yield  » 3 T with  no  armature 
reaction. 

Similar  effects  due  to  armature  reaction  occur  for  rectangular  geometries. 
However,  It  Is  easy  to  comoensate  rectangular  generators  by  returning  the  load 
current  through  the  magnet  air  gap.  For  an  annular  geometry,  compensation 
would  mean  running  the  load  current  through  the  fluid  gap  between  the  electrodes, 
and  this  would  eliminate  one  prime  advantage  of  the  annular  geometry— the  lack 
of  Insulating  walls  between  the  two  electrodes. 
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Table  IV. 6 

Powers  and  Efficiencies  With  and  Without  Armature  Reaction, 
r^  = 1.0  m,  Rq  =»  1.1,  L = 1 m,  uz  » 15  m/s,  Vi  ■ 3.903, 

F = 0.91,  a * 2.62  x 10^  mhos/m,  and  R^,  - 4.939 


Armature 

Reaction 

Ignored 

Connection 

Downstream 

Upstream 

Both  Ends,  L'  = 0.8  m 

II  [MA] 

4.445 

* 124.7 

0.8936 

2.399 

Pg  [MW] 

17.92 

502.8 

3.602 

9.673 

P„  [MW] 

19.17 

2,759 

3.824 

10.61 

[MW] 

1.14 

2,256 

0.1158 

0.835 

P/[MW/m''’] 

27.16 

762.1 

5.46 

14.66 

Ap/L  [N/m^] 

1.937  X 10® 

’ 2.789  X 10® 

3.865  X 10® 

1.073  X 10® 

Hg  [%] 

93.5 

18.2 

94.2 

91.1 

(z-0)  [T] 

3.0 

3.0 

2.821 

2.823 

(z-L)  [T] 

3.0 

27.94 

3.00 

3.303 

*Pj  * useful  power  density 


'j 

•1 

ii 

i 


I 
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4 


Bi  (T) 


with  *rw«ture  Reaction 


V.  CONCLUSIONS  AND  FUTURE  WORK 


The  ambient-temperature  NaK-nItrogen  facility  has  been  successfully  revised 
to  handle  higher  velocities  and  pressures,  and  satisfactory  operation  has  been 
obtained.  The  Improved  control  features  make  the  system  easier  to  use.  New 
experiments  with  the  second  dlverglng-'hannel  generator  showed  that  the  power 
density  was  higher,  and  that  there  appeared  to  be  e tendency  for  the  efficiency 
to  Improve  more  with  Increasing  quality  at  higher  velocities. 

The  program  to  produce  and  evaluate  liquid-metal  foams  was  Initiated. 

Equipment  was  developed,  screening  tests  made  with  mercury  while  waiting  for 
the  high-purity  NaK,  and  a nlan  established.  The  future  program  Is  described 
briefly  below. 

The  evaluation  of  an  annular  generator  geometry  was  completed.  Although 
this  geometry  appears  advantageous  because  there  Is  no  need  for  Insulating  walls, 
the  cylindrical  (electrode)  shells  are  Ideal  pressure  vessels,  and  the  magnetic 
field  Is  contained  entirely  within  the  generator;  there  are  also  serious  limitations. 
The  limitations  are: 

1.  The  geometry  Is  Inherently  suited  to  low  voltage  high-current  operation, 
not  to  the  higher  terminal  voltages  needed  as  Input  to  acycl 1c  motors  or  to 
efficient  and  cheap  dc  to  ac  converters. 

2.  Flow  reversal  and  Increased  viscous  loss  can  occur  at  the  power  densities 
desired  for  power  systems. 

3.  Armature  reaction  effects  are  significant  and  compensation  methods 
appear  Impractical. 

Although  these  limitations  were  known  before  the  analysis  described  In  Section  IV 
was  performed,  the  analysis  was  necessary  to  quantify  the  effects.  Based  on  this 
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analysis,  this  annular  geometry  does  not  appear 
pump  operation. 

The  desired  next  steps  In  the  liquid-metal 
below. 

V.l  Higher  Velocities  and  Pressures 

The  motivations  for  tests  at  higher  liquid 
slip  loss  due  to  the  relative  ohase  velocities, 
demonstrated  that  this  loss  can  be  decreased  by 
and  to  test  the  generator  at  velocities  typical 
the  experiments  to  date  exceot  for  a few  recent 
been  limited  to  about  12  m/s.  The  revised  facility  allows  generator  operation 
at  higher  liquid  velocities  or  flow  rates  and  higher  pressures.  Testing  of  the 
existing  diverging  channel  will  continue  at  liquid  flow  rates  up  to  about 
10.9  kg/s  (200  gpm),  versus  the  previous  limit  of  about  6 kg/s  (110  gpm).  This 
generator  was  designed  for  lower  velocities  and  has  an  operating  pressure  range 
up  to  about  0.72  MPa  absolute  (90  psig),  so  that  operation  at  higher  velocities 
will  require  a reduced  magnetic  field  strength. 

A new  generator  will  be  designed  and  built  to  utilize  the  full  loop 
capability  of  10.9  kg/s  and  1.48  MPa  absolute  (200  psig).  This  will  allow 
high  velocities  with  high  field  strengths,  and  variation  of  the  generator 
pressure  range  with  a maximum  generator  exit  pressure  of  1.14  MPa  absolute 
(150  psig).  This  generator  will  have  a different  taper  to  fit  the  different 
design  pressure  ratio.  Also,  the  maximum  void  fraction  In  the  generator  should 
be  Increased— this  will  be  aided  by  a reduced  velocity  slip  ratio. 


attractive  for  either  generator  or 

MHO  program  are  described 


velocities  are  to  reduce  the 
as  experimental  research  has 

3 

Increasing  the  liquid  velocity, 
of  practical  generators.  In 
runs,  the  liquid  velocity  has 


1 

I 


V.2  Foams 

The  Initial  steps  of  the  program  to  evaluate  foams  In  MHD  generators  were 
described  In  Section  I I. 3.  In  the  next  year  attention  will  be  focussed  on 
Identifying  surface-active  agents  suitable  for  NaK,  and  plans  should  be  developed 
for  testing  NaK  foams  In  a generator.  Bubbling  tests  with  NaK  and  various 
additives  will  be  part  of  the  screening  program,  along  with  measurements  of 
surface  tension,  viscosity,  and  wetting  oropertles.  NaK  was  chosen  for  the 
Initial  tests  because  of  extensive  experience  with  NaK  In  generators,  and  thus 
the  ability  to  see  the  change  In  generator  performance  due  to  the  use  of  a 
foam  In  minimum  time  and  for  minimal  cost.  If  the  NaK  test  results  are  positive, 
then  effort  will  be  directed  toward  finding  foams  suitable  for  use  In  practical 
MHD  systems.  Methods  of  separating  the  gas  and  liquid  phases  after  the  MHD 
generator  for  a foam  must  also  be  considered. 

V.3  Improved  Local  Diagnostics 

Improved  local  diagnostics  are  needed  to  obtain  a better  understanding  of 
two-phase  flows  In  a magnetic  field.  The  data  to  be  measured  Include  direct 
local  values  of  liquid  and  gas  velocities,  void  fractions,  bubble  diameters, 
pressures,  and  temperatures.  Measurement  techniques  to  be  evaluated  Include: 

(1)  conductivity  orobes,  both  single  and  multiple,  for  void  fractions,  bubble 
sizes,  and  bubble  velocities;  (2)  hot-film  probes  for  liquid  and  gas  velocities, 
void  fractions,  and  bubble  sizes;  and  (3)  microthermocouples  for  gas  and  liquid 
temperatures.  These  techniques  have  been  used  successfully  for  air-water  two-phase 
flows,  but  have  not  been  tried  with  liquid  metals  or  magnetic  fields.  Special 
electronics  and  data-processing  techniques  may  be  required. 
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This  local  data  Is  needed  to  develop  and  evaluate  Improved  analytical  models 


for  two-phase  flows,  Section  V.4.  Thus,  a program  to  develop  the  measurement 
techniques  should  be  Initiated.  First  tests  should  use  the  air-water  facility 
described  In  Reference  4 before  Incorporation  Into  the  ambient-temperature 
NaK-nItrogen  loop. 

V.4  Analytical  Models  for  the  Two-Phase  Flow  In  the  Generator 

Further  analytical  studies  on  the  two-ohase  flow  In  the  MHD  generator  should 
be  carried  out  In  parallel  with  the  above  efforts.  The  experimental  data, 
especially  that  from  the  Improved  local  diagnostics  as  It  becomes  available, 
will  be  used  In  developing  and  evaluating  these  models.  The  primary  goals  are 
a better  understanding  of  and  model  for  the  two-phase  flow  at  high  void  fractions 
In  a practical  generator.  The  conditions  for  transitions  between  bubbly, 
churn-turbulent,  slug,  and  annular  flows  at  high  void  fractions  may  be  explored. 

Of  particular  Importance  are  the  effects  on  generator  performance  of  velocity, 
void  fraction,  and  conductivity  distributions  across  the  channel,  I.e.,  parallel 
to  either  the  applied  magnetic  field  or  the  current.  These  have  not  been 

3 

Included  In  any  model  to  date  except  for  a simple  boundary-layer  analysis. 

i 

V.5  Imperfect  Compensation 

j 

Most  analyses  of  two-phase  MHD  generators  have  assumed  perfect  compensation,  ! 

I.e.,  the  generator  current  Is  returned  through  the  air  gap  In  the  magnetic  ] 

i 

circuit  In  such  a fashion  that  the  applied  magnetic  field  Is  not  changed  by  the  j 

1 

generator  current.  This  may  not  be  true  In  practice  because  of  end  currents,  | 

non-uniform  current  distributions  In  the  generator  and/or  compensating  conductors,  j 

1 

J 

I 
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and  compensating  conductors  that  are  longer  than  the  electrodes  as  in  the 
present  generator.  The  development  of  an  analytical  model  for  the  generator 
with  imperfect  compensation  has  been  initiated.  Future  plans  include  experi 
mental  verification  of  the  model. 
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APPENDIX  A 

SUMMARY  OF  EXPERIMENTAL  GENERATOR  DATA 
AUGUST  1975  TO  SEPTEMBER  1976 


Date 

Run 

Number 

Mass  Flow  Rate 

NaK  ^*2 

• * 

[kg/s]  [kg?s] 

Magnetic 

Flux 

Density 

B 

[T] 

Load 

Resistance 

•^L 

[mft] 

Load 

Voltage 

[Ji 

Pressure 
Inlet  Exit 

Pin  ^out 
[psia]  [psia] 

Power 

Output 

% 

[kW] 

Average 

Void 

Fraction 

8-12-76 

1 

6.025 

0 

1.245 

00 

.509 

63.3 

37.5 

0 

0 

2 

6.115 

.0125 

1.246 

00 

.717 

57.7 

34.2 

0 

.277 

3 

6.025 

.03175 

1.244 

- 

.932 

57.1 

35.2 

0 

.480 

4 

6.025 

.0507 

1.245 

00 

1.102 

59.5 

36.8 

0 

.558 

5 

6.070 

.0655 

1.245 

00 

1.257 

60.7 

36.0 

0 

.652 

8-16-76 

1 

12.140 

0 

1.245 

00 

.996 

88.7 

31.0 

0 

0 

2 

12.095 

.0135 

1.235 

- 

1.316 

74.3 

26.4 

0 

.184 

3 

11.560 

.0359 

1.234 

00 

1.784 

78.4 

31.0 

0 

.450 

8-18-76 

1 

6.160 

.0351 

1.235 

.2706 

.867 

88.4 

34.0 

2.78 

.495 

2 

6.003 

.0516 

1.238 

.271 

.932 

92.6 

36.3 

3.21 

.587 

3 

6.025 

.0672 

1.235 

.271 

1.010 

93.5 

33.7 

3.76 

.678 

4 

6.025 

.0130 

1.240 

.2703 

.665 

80.3 

31.5 

1.63 

.273 

5 

6.025 

0 

1.253 

.264 

.470 

81.1 

35.0 

.837 

0 

6 

6.025 

.0654 

1.235 

.270 

.931 

77.0 

18.5 

3.21 

.672 

3-19-76 

1 

6.038 

0 

1.248 

.2644 

.476 

65.0 

18.3 

.856 

0 

2 

6.025 

.01396 

1.233 

.2709 

.852 

70.7 

15.5 

2.70 

.524 

3 

6.025 

.0356 

1.230 

.2704 

.725 

65.4 

16.4 

1.94 

.388 

4 

6.025 

0 

1.25 

00 

.492 

45.1 

18.1 

0 

0 

5 

11.825 

0 

1.265 

00 

.982 

64.3 

15.5 

0 

0 

6 

6.070 

.01231 

1.236 

00 

.827 

39.2 

17.5 

0 

.435 

7 

6.025 

.03078 

1.235 

00 

1.102 

36.8 

21.7 

0 

.608 
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Date 

Run 

Number 

Mass  Flow  Rate 
NaK  J|2 

[lcg?s]  [kg?s] 

Magnetic 

Flux 

Density 

B 

[T] 

Load 

Resistance 

[mh] 

Load 

Voltage 

[Hi 

Pressure 
Inlet  Exit 

**10  ^out 

[psial  [psia] 

Power 

Output 

Pe 

[kW] 

Average 

Void 

Fraction 

BT 

8-24-76 

1 

6.02 

0 

1.223 

.240 

.466 

81.6 

35.4 

.903 

0 

2 

5.98 

.06838 

1.222 

.2468 

.937 

93.0 

26.6 

3.56 

.621 

3 

6.07 

.05049 

1.227 

.2475 

.885 

91.5 

35.3 

3.13 

.547 

4 

6.025  .03512 

1.227 

.2471 

.824 

87.4 

33.7 

2.75 

.477 

8-25-76 

1 

6.025  .07792 

1.225 

.2461 

1.003 

95.3 

37.6 

4.08 

.542 

2 

12.05 

.06697 

1.26 

Oft 

1.988 

83.2 

34.0 

0 

.394 

3 

6.025  .01372 

1.225 

.2441 

.662 

87.5 

36.0 

1.80 

.161 

4 

6.07 

0 

0 

«0 

0 

31.8 

31.7 

0 

0 

9-2-76 

1 

12.09 

0 

0 

«» 

0 

32.0 

37.5 

0 

0 

2 

6.025  .03564 

.805 

OB 

.713 

42.7 

35.1 

0 

.574 

3 

5.98 

.05139 

.807 

e» 

.822 

45.0 

35.3 

0 

,637 

4 

6.07 

.01393 

.806 

€» 

.544 

43.7 

34.3 

0 

.390 

5 

6.025 

0 

.808 

- 

.313 

44.3 

32.2 

0 

0 

6 

12.09 

.05136 

.805 

00 

1.439 

56.3 

30.8 

0 

.438 

9-9-76 

1 

12.09 

0 

0 

0 

33.3 

38.0 

0 

0 

2 

5.98 

.05797 

.801 

- 

.866 

40.5 

36.4 

0 

.734 

9-13-76 

1 

18.03 

0 

0 

«» 

0 

44.8 

47.7 

0 

0 

2 

5.98 

.09896 

1.35 

.261 

1.008 

97.7 

34.8 

3.89 

.581 

9-15-76 

1 

12.14 

0 

.805 

.251 

.617 

78.0 

36.3 

1.52 

0 

2 

11.91 

.09533 

.6 

.260 

1.160 

74.5 

28.4 

5.18 

.601 

3 

11.96 

.07093 

.595 

.259 

1.018 

71.7 

33.2 

4.00 

.511 

4 

12.095  .03775 

.603 

.259 

.825 

65.5 

30.3 

2.63 

.389 

5 

12.21 

0 

.620 

.251 

.472 

54.2 

28.7 

.887 

0 

6 

6.025  .04863 

.605 

.259 

.511 

51.0 

34.3 

1.01 

.546 
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Date 

Run 

Number 

Mass  Flow  Rate 
NaK  J|2 

[kg^s]  [kg?s] 

Magnetic 

Flux 

Density 

B 

[T] 

Load 

Resistance 

h 

[mn] 

Load 

Voltage 

Pressure 
Inlet  Exit 

Pin  Pout 

[psia]  [psia] 

Power 

Output 

Pe 

[kW] 

Average 

Void 

Fraction 

S’ 

9-17-76 

1 

6.151 

0 

.6 

.250 

.232 

47.0 

36.1 

.215 

0 

2 

6.002 

.01261 

.6 

.255 

.344 

49.5 

36.4 

.464 

.317 

3 

5.98 

.03318 

,595 

.258 

.451 

52.0 

36.8 

.788 

.520 

4 

6.07 

.06330 

.595 

.259 

.539 

53.8 

35.1 

1.12 

.610 

5 

11,915 

.07023 

.588 

.373 

1.068 

64.2 

32.0 

3.06 

.541 

9-20-76 

1 

5.980 

0 

.599 

.238 

43.9 

36.6 

0. 

0 

2 

6.025 

.0135 

.613 

00 

.403 

41.8 

35.4 

0 

.359 

3 

6.025 

.0348 

.598 

oo 

.543 

41.7 

35.8 

0 

.585 

4 

6.025 

.0502 

.598 

00 

.638 

42.9 

34.4 

0 

.637 

5 

5.98 

.06378 

.602 

00 

.695 

43.7 

35.5 

0 

.726 
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